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Submarine clay sediments exhibit time dependent stress-defor- 


mation characteristics and can be described by a viscoelastic 


relaxation modulus. A method is developed for determination of 


the relaxation modulus in shear which utilizes viscoelastic theory 


and torque versus rotation data from vane shear tests. The modu- 


lus, G(t), is described by the power law: 

g(t) = Gt 
where t is time, and G, and n are constants. 4G, and n are deter- 
mined from vane tests on deep sediment core samples. G, is rota- 
tion angle dependent. The validity of the procedure is supported 
by predictions of in situ vane test torque versus rotation curves 
which agree favorably with data obtained using an in situ vane 


device. G, correlates with maximum vane shear strength and n 


correlates with liquidity index. Additional work is indicated. 
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INTRODUCTION 


In recent years foundation engineers have realized the impor- 
tance of developing a reliable and accurate means of characterizing 
the stress-deformation behavior of submarine sediments. The only 
parameters available for design use have been the so-called engi- 
neering properties of the sediment, as determined from laboratory 
tests. These properties include ultimate shear strength (cohesion), 
Atterberg limits and other measures upon which much of the practice 
of soil mechanics is based. These parameters yield no quantitative 
information as to the stress-deformation behavior of the sediments 
prior to failure. In particular, they provide no means of 
characterizing time (and frequency) dependence of the sediment 
shearing resistance. Time dependence of the stress-strain relation- 
ship is particularly important in a large class of problems where 
dynamic, cyclic or vibration loadings exist, or in long term 
(creep) deformation analysis. 

Sea floor instability has been recognized as a serious engi- 
neering problem (38)*. Review of previous efforts at prediction 
and analysis of instability of deltaic sediments will illustrate 


the need for characterizing the time dependency of the sediment. 


*Numbers in parenthesis refer to references in Appendix I. 
(The citations follow the style of the Journal of Soil Mechanics 


and Foundations Division, American Society of Civil Engineers. ) 





Slides are known to have occurred in areas off the Mississippi 
River Delta of very slight slope, where conventional statical 
equilibrium analysis (using ultimate Strength) yields factors 

of safety well above 1.0. It has thus been concluded that these 
Slides are not entirely gravity controlled but are triggered by 
some mechanism causing the sediment to become unstable. Henkel 
(20) first produced a simplified model to indicate that bottom 
instability could be induced by cyclic hydrostatic bottom pres- 
sure variations caused by storm waves. Subsequent theoretical 
and experimental investigations have tended to verify that storm 
waves can induce instability and more complex methods of analysis 
have been proposed (42, 26). Wright and Dunham (42) developed 

a finite element model for the state of stress and deformation 
in the sediment when subjected to wave induced cyclic shear 
stresses. They treated the problem as one of static elasticity 
and assumed a hyperbolic stress-deformation relation. This, 
together with measured values of ultimate shear strength and an 
assumed value of Poisson's ratio, allowed estimation of an elastic 
shear modulus (non time-dependent). Gravity is neglected in the 
model and no information regarding the actual failure mechanism 
is obtainable. Mitchell, Tsui and Sangrey (26) set forth a 
method of analysis based on a gravity controlled failure which 
occurs because of wave induced reduction in soil shear strength. 


As with the finite element analysis, the Strength reduction 





approach is dependent on a measure of the ultimate shear strength 
of the sediment. No stress-strain characteristics are used and 
no information concerning the actual mechanism of failure is 
obtained. Additional experimental and theoretical treatments of 
submarine slides were recently reported (8, 1). These serve to 
further confirm that storm waves are the likely trigger mechanism 
but have added nothing toward an understanding of the mechanism 
of failure. 

To fully understand the phenomenon a broad approach must be 
taken. The following elements are deemed essential to a realistic 


analysis: 


1. Prediction of the state of the sea surface and the result- 


ing particle motions and the dynamic pressures in the 
water column. 

2. Determination of the dynamic stress-deformation charac- 
teristics of the bottom and sub-bottom sediments. 

3. Analytical determination of the dynamic response of the 
bottom and sub-bottom material prior to failure and the 
interaction of this response with the water motion. 

4. Determination of the sea state (wave conditions) that 
will produce complete bottom instability, and prediction 
of the subsequent characteristics of the failed mass. 

Available analytical and statistical methods (21) can provide 


a reasonably accurate solution to element 1 above. Element 3 has 





been partially treated by Gade (14) although this work dealt more 
with the effect of the bottom on the waves rather than the waves 
on the bottom. The works previously mentioned (1, 8, 26, 42) 
have attempted to provide answers relative to element 4; however, 
without sediment stress-deformation characteristics under dynamic 
loading and a realistic model of the combined motion of the water- 
sediment system (elements 2 and 3 above), these analyses have not 
provided a general slope instability theory. Schapery (34) has 
proposed a theoretical wave analysis of the combined motion of 
the water-sediment system which may provide this realistic model 
and lead to characterization of the onset and occurrence of 
instability. This work is currently in progress. Determination 
of the dynamic stress-deformation characteristics of the in situ 
sediments (element 2 above) remains. 

With the foregoing in mind the general objective of this 
study was to investigate the in situ dynamic stress-deformation 
characteristics of sea floor clay sediments and to determine a 


means of quantifying these characteristics. The word in situ is 





a key word here. It is well known that sediment samples recovered 
from the sea floor undergo considerable and unpredictable distur- 
bance due to mechanical sampling disturbance, expansion of gases 
in the sample, and other factors (28). Even with the most careful 
sampling procedures, the material properties obtained from labora- 


tory tests of these samples cannot be called true measures of 


these properties for the in situ material. The implication is, 
then, that determination of in situ material relationships should 
be based upon an in situ test conducted, insofar as possible, in 
the natural sediment environment. 

The only known in situ soil strength measuring device capable 
of obtaining in situ sediment shearing resistance data at deep 
penetrations in bottom sediments is the wire line down hole vane 
shear device in use by McClelland Engineers, Inc. (9). This 
device has been used extensively for ultimate shear strength 
measurements. In addition, the vane in this device can be rotated 
at different rates and a continuous readout of torque (measure of 
sediment resistance to shear) versus rotation angle (measure of 
deformation) is available. 

At the outset then, it was decided to pursue development of a 
procedure to utilize a down hole vane shear device to quantify a 
dynamic stress-deformation relation for sea floor clays. Schapery 
(34) proposed a method of obtaining such a relation from torque 
versus rotation vane shear data which treats the sediment as a 
nonlinear viscoelastic material. Specifically the proposed pro- 
cedure results in the computation of a time dependent (viscoelas- 
tic) relaxation modulus in shear which can be generalized into a 
complex shear modulus. 

This research, then, was undertaken with the following 


specific objectives: 





To investigate (by literature review) the factors affect- 
ing the stress-deformation characteristics of sea floor 
clays and to review previous studies of the use of vane 
Shear apparatus. 

Determine the feasibility of characterizing sea floor 
clays by a viscoelastic shear modulus and of computing 
the modulus from vane shear data. 

To investigate the use of the down hole vane shear device 
to obtain the required data for shear modulus determina- 


tion. 





BACKGROUND 


No published information relating vane shear test data to 
viscoelastic sediment properties has been found in the litera- 
ture. Several recent studies have characterized sea floor 
sediments as viscoelastic materials but few have quantified a 
relationship; although viscoelastic theory and experimental pro- 
cedures are highly developed. The vane shear test has been widely 
studied as a means to determine the ultimate shear strength of 
clays and recently has been extended to in situ use for sea floor 
clays and more sophisticated laboratory studies. 

In this section the basic nature and mathematical treatment 
of viscoelastic materials will be briefly reviewed; the literature 
characterizing clays as viscoelastic materials will be reviewed; 
and the proposed physical basis for the viscoelastic stress-defor-. 
mation characteristics of sea floor sediments will be summarized. 
In addition, a review of the use of the vane shear test will be 
made. 

Characterization of Viscoelastic Materials. - A material is 
said to be viscoelastic if stress in it is a function of current 
and/or past values of strain and strain rate. In more general 
terms, stress is a function of strain history. Ferry (12), Fung 
(13), Pipkin (30), among others, have presented detailed treat- 


ments of the fundamentals of viscoelasticity theory. In the 


theory of linear viscoelasticity the material is often viewed as 
possessing properties of both linearly elastic (Hookean) solids 
and Newtonian fluids, and they are represented by combinations 
of the familiar spring and dashpot models (30, 44). Mathemati- 
cally each Hookean element deforms in simple shear stress condi- 
tions according to: 

CREO VY (G) ene. ee sc ese kw ee (1) 
where y(t) is the shear strain, t(t) is the corresponding stress 
and G is the elastic modulus in shear. Each Newtonian element 
deforms according to: 

26) 6 TEI... ee eee (2) 
where n is the fluid viscosity and y(t) is the strain rate. In 
addition to the two fundamental model elements for representing 
viscoelastic materials, some investigators have employed the non- 
linear yield stress element, which is represented by a block, 
resting on a surface, which requires some minimum (yield) stress, 
Ty, to set it in motion. Fig. |] depicts the elementary model 
elements along with several common rheological models which have 
been employed for describing soil behavior. The differential 
equations of the associated stress-strain relations in shear are 
also shown. Yong and Warkentin (44) present additional models 
and develop the associated stress-strain equations. Suklje (39) 
presents an extensive treatment of the rheological nature of soils, 


including development of viscoelastic theories. 
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FIG. 1. - SIMPLE RHEOLOGICAL MODELS 





These simple models are useful in visualizing material] behav- 


ior, but, for most real materials, are too simple to provide 
accurate predictions of materia] behavior unless a large number 
of elements are used. Therefore, a more general mathematical 
statement for viscoelastic stress-strain behavior is desired. 
Linear viscoelastic materials, by definition, obey the law of 
superposition which says the response due to two simultaneous 
inputs is equal to the response due to one input alone, plus 
‘the response due to the other. These materials also obey the 
proportionality relation, that, at any t, the response due to 
an input times a constant equals the constant times the response 
due to input alone. For these materials the stress (or strain) 
response due to any continuous time dependent strain (or stress) 
input can be obtained by summing the responses due to an arbi- 
trarily large number of input "steps" from the time of the 
beginning of load application (for nonaging materials) to any 
time of interest (30). For the case of a series of shear strain 
input steps we thus may obtain: 
n 

T(t) : eal i, ayy eee eee (3) 
where: G(t) is the shear modulus, a functional of time, n is the 
number of strain steps, AY, is the strain change at step a and 
t - t is the time elapsed since application of strain step «. 


By passing this sum to the limit we obtain the so called heredity 





or convolution integral for a linear, nonaging viscoelastic mate- 
rial undergoing any (time dependent) simple shearing strain input 


beginning at the time t = 0 (30). 


where G(t) is the time dependent modulus and t' (taking the place 
of t. Eq. 3) is a dummy variable of integration. The convolution 
equation, (Eq. 4), based on the linearity assumption, forms the 
basis for much of the mathematical treatment of viscoelasticity, 
and can be applied to the solution for stress or strain responses, 
given the function G(t), or as a means of finding the response 
function, G(t), from experimental stress-strain data. Eq. 4 will 
be utilized in later developments in this study. 


Viscoelasticity in Sea Floor Sediments - Physical Basis. - 
Investigations by Singh and Mitchell (37) and Mitchel et al. (25), 


among others, have provided experimental evidence that creep defor- 


mation in clay soils can be treated by viscoelastic stress-strain- 
time functions. On the basis of these results Schapery (32, 33) 
has included clays among such diverse materials as polymers, 
metals and biological tissue in characterizing nonlinear visco- 
elastic materials. There are many additional examples in the 
literature of the treatment of clays, in various modes of defor- 
mation, using viscoelastic rheological models or analyses (29, 


35, 39, 43). In addition, several recent investigations have 
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specifically characterized sea floor clays as viscoelastic mate- 
rials (36, 27, 6). Sherif et al. (36) and Monney (27) proposed 
the nonlinear model shown in Fig. 2 and showed creep test data 

on submarine sediments which were in close agreement with those 


predicted by the model. Monney (27) has pointed out that this 


G5(e) 
Go) To 
T) (o) — 
20 O, 


FIG. 2. - NONLINEAR RHEOLOGICAL MODEL FOR SEA FLOOR CLAYS 
PROPOSED BY SHERIF ET AL. [FROM SHERIF ET AL. (35)] 


model is consistent with the work of Cohen (6) which indicated 
that nondispersed clay sediments behaved as nonlinear viscoelastic 
materials. Other experimental evidence supporting the treatment 
of submarine clays as viscoelastic materials is presented by 
Carpenter et al. (5) who used shear viscometer data on undis- 
turbed samples of Mississippi Delta clays to obtain a constitu- 


tive equation for stress which involved strain and strain rates. 





Ue, 


The physical mechanism of the stress-deformation response of 
clays is not well understood. However, a review of the probable 
causes of rate dependent behavior in submarine sediments is appro- 
priate. Scott (35) and Wu (43) present general reviews of the 
factors affecting clay structure and stiffness. According to 
Sherif et al. (36) viscoelasticity in submarine sediments is 
largely a result of very high water contents and the type of clay 
structure. Salt water deposited clays can be thought of as com- 
posed of mineral platelets arranged in the so called flocculent 
or card house structure, such as depicted in Fig. 3, although 
several different proposals have been advanced for the exact 


structure (24). It is believed that, during the marine 





Flocculent Nispersed 


FIG. 3. - DEPICTION OF FLOCCULENT AND DISPERSED CLAY STRUCTURES 
[FROM WU (43) ] 


depositional process, the individual clay particles acquire strong 
positive charges along their edges due to high hydrogen ion concen- 
trations in these regions (43). In addition, due to isomorphous 


substitution (43) the flat sides of the platelets take on negative 








surface charges. The result is that, as these particles settle, 
there are strong edge to face attractive forces which overcome 
weaker repulsive forces and cause the particles to come into con- 
tact or bond together in the edge to face or card house or floccu- 
lent structure (43). Lambe (22) has characterized the attractive 
forces as largely van der Waal‘s forces. 

The water in the interstital spaces exists in two states, 
pore water and absorbed water. The pore water is "free" in the 
interstital spaces while the adsorbed water is a layer or film 
of water molecules held tightly to each clay particle by hydro- 
gen bonding and other forces. According to Grim (16) the 
adsorbed water is actually more dense, and much more viscous 
than the pore water. It must be noted that the above descrip- 
tion of clay structure and formation is highly simplified. Com- 
plex ionic interactions are involved in the depositional process 
which depend on the type of clay mineral involved and the ionic 
content of the water (16); however, treatment of these in detail 
is beyond the intent of this discussion. 

From this simplified model of clay structure the process of 
deformation under shear can be considered and some of the factors 
affecting the viscoelastic nature of submarine clays can be 
identified. When the flocculent clay structure is subjected to 
Shear the bonds between the particles are first strained. As 


the bonding forces are progressively overcome along a shear 
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surface, and shearing continues, the viscous shear resistance of 
the pore water becomes important as the clay particles along the 
shear surface tend to orient themselves into a more dispersed or 
oriented structure (see Fig. 3, page 13) (35). 

Interparticle contacts along the shear surface are broken 
as the particles move to the dispersed structure. This may have 
several effects. Part of the normal stress across the failure 
surface will be shifted to the pore water causing an increase in 
pore water pressure. The increased pore water pressure tends to 
force water out of the strained area. If sheared rapidly the 
water cannot escape (undrained condition) and the result is 
decreased shear resistance. If sheared very slowly (drained 
condition), so that the pore water can escape the effect of the 
shifting of the particles to a dispersed structure will be less. 


The shear resistance will then tend to be greater. The rate of 


migration (drainage) of the pore water in the clay is, of course, 


a function of its viscosity. 
Leonards and Girault (23) have presented evidence that the 
adsorbed water may also play a role in the viscoelastic behavior 


of clays. 


Summarizing, viscoelastic behavior in sea floor clays may be 


attributed to a combination of clay structural changes, shearing 
and migration of the pore water, and effects of the adsorbed 


water. These phenomena are controlled by the type of clay 
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mineral, the ionic constituents of the water at the time the clay 
was deposited, the void ratio (i.e. relative amount of water in 
the pores for saturated sediments), the viscosity of the pore 
water, and the rate of shearing. 

Experience With the Vane Shear Test. - Since the classic 
work of Cadling and Odenstad (4) the vane shear test for shear 
strength of clays has become a widely used test both in situ 
and in the laboratory. The conventional vane test consists of 
turning a four bladed (usually rectangular blades) vane in the 
soil and measuring the torque required to shear the soil around 
the cylinder of revolution of the vane. The basic test is 


depicted in Fig. 4. Until recently almost all vane shear 






FOUR BLADED 
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FIG. 4. - BASIC VANE SHEAR TEST 





instruments were hand operated and the torque was measured by 
noting the deflection of a calibrated spring integral to the 
vane shaft. This discussion is intended to briefly summarize 
the theory and use of the vane shear test and to describe recent 
investigations of its use on submarine sediments, both in situ 
and in the laboratory. 

Cadling and Odenstad (4) presented the most complete mathe- 
matical treatment of the vane test to date and provided the first 
substantial experimental evidence of its applicability in clays. 
Their experiments provided evidence that: 

1. The failure surface caused by rotating a four bladed 
vane in clay is closely approximated by a right 
circular cylinder. 

2. The effect of progressive failure (failure starting 
in front of the edge of the blades and propogating 
along the cylindrical surface) was negligible in 
ultimate strength tests. 

3. The strength of the clay was affected by the vane 
rotation rate. 

4. A vane length to diameter ratio of about two was most 
suitable for in situ tests. 

Gibbs (15) offered further evidence that the failure surface 

is cylindrical. Eden and Hamilton (10) and Gibbs (15) showed 


that vane insertion causes minimal distrubance effects. 





The calculation of the shear strength from vane tests using 


rectangular blades is normally accomplished using the equation: 


i C, ( 5 Bo (5) 
where 
Lee = maximum vane shaft torque measured to failure 
Ch = ultimate soil shear strength (also called cohesion) 
d = vane diameter 
L = vane length 


Eq. 5 is obtained by summing the torque due to the shear 
resistance of the clay over the cylindrical failure surface, 
and the top and bottom. 

The vane test is particularly applicable to very soft 
clays for which it is difficult to obtain and test undisturbed 
samples by other soil strength tests (4, 41). L. Casagrande 
first reported use of the in situ vane shear test in underwater 
sediments (4). He used the test in river sediments at water 
contents above the liquid limit. These and later tests (4) 
indicated the vane shear determined strengths were 1-2 times 
the strengths determined on laboratory samples using other tests. 
Gray (17) compared strengths from field vane tests to those 
obtained from laboratory tests of samples obtained in thin 
walled core tubes. The clay was sensitive and low strength. 


The laboratory vane shear strengths were less than the in situ 
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values but were greater than the strengths as determined from 
other laboratory tests. The lower laboratory values were attri- 
buted to sampling disturbance (17). 

Fenske (11) conducted the first comprehensive study of the 
use of the vane test on submarine sediments. His tests were 
conducted in Mississippi River delta sediments off the coast of 
Louisiana. In situ vane shear strengths were measured to depths 
of 24] feet where the water depth was 66 feet. The device used 
was a combination vane shear device and borehole drilling head 
as shown in Fig. 5. Fenske conducted extensive comparisons of 
in situ vane, minature laboratory vane from undisturbed cores and 
unconfined compression strengths. Fig. 6 summarizes the compari- 
son and indicates that the act of sampling and raising submarine 
sediments to the surface changes the strength properties. Hall 
(18) demonstrated the use of the more refined Swedish Foil 
Sampler to reduce sample disturbance, but this device is complex 
and is not easily applicable to ocean work. The occurrence of 
disturbance due to submarine sampling and recovery is largely a 
result of mechanical sampling disturbance and the expansion of 
dissolved gases as the hydrostatic pressure is reduced (28). 

With the increasing requirement to design foundations for 
sea floor structures several additional recent investigations of 
in situ vane testing have been conducted (31, 7, 9). Richards 


et al. (31) and Demars and Taylor (7) were concerned with 
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FIG. 5. - DOWN HOLE SUBMARINE VANE TEST AND BORER OF FENSKE 
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FIG. 6. - COMPARISON OF IN SITU AND LABORATORY ULTIMATE SHEAR 
STRENGTHS [AFTER FENSKE (11)] 


evaluation of bottom resting devices capable of obtaining vane 
Shear strengths to sediment depths of only a few feet. These 
Studies are further evidence of the applicability of the vane 
test and tend to confirm the comparisons obtained by Fenske. 
Doyle et al. (9) have demonstrated a remotely controlled wire 
line vane test device, developed by McClelland Engineers, Inc., 
for down hole vane measurements in undersea clays. The device 

is lowered inside the borehold casing and can be used alternately 


with the 2.25 inch wire line core sampler used extensively 





offshore. This allows practica!iy continuous in situ vane testing 
and "undisturbed" sampling for she full length of a bore hole. 


The device and sequence of operation is depicted in Fig. 7. It 
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FIG. 7. - DOWN HOLE WIRE LINE SUBMARINE VANE SHEAR DEVICE OF 
McCLELLAND ENGINEERS [FROM DOYLE ET AL. (9)] 


can be operated from a floating platform and is considered by the 
author to be the only practicable deep penetration sea floor test 
device presently available. The wire line vane device is capable 
of several different vane rotation rates and of measuring torque 


versus rotation angle data. (This fact will be useful later in 


this study. ) 
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Coincidental with the development of in situ vane tests there 
has been a recent trend toward sophistication of the laboratory 
vane instruments for ag floor sediment sample tests. Precision 
and sensitivity beyond the capabilities of hand cranked, spring 
loaded torque measurement devices is required in the study of 
sea floor sediments. In addition, the recognition of the impor- 
tance of the rheologic nature of sea floor clays has prompted the 
development of vane shear instruments which provide a constant 
readout of torque versus rotation (27). The most recent instru- 
ments employ controlled speed motorized vanes and measure torque 
using strain gauged torque transducers (27, 19). These devices 
avoid the problems of inaccurate determination of actual vane 
rotation angle and speed, and the build up of strain energy in 
the load measuring spring of the conventional devices is elimi- 
nated. Furthermore, a continuous record of torque versus vane 
rotation is obtainable. It shall be later demonstrated that 
these features are essential to the vane shear determination of 


the viscoelastic properties of submarine clays. 
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THEORETICAL BASIS FOR EXPERIMENT 


One objective of this research was to develop a procedure to 
determine a viscoelastic shear modulus for sea floor sediments. 
Such a procedure can be only as good as the samples tested. As 
noted in the previous section the mere act of sampling and recov- 
ering samples seriously disturbs them. Thus a determination based 
upon in situ measurements would be most desirable. As a matter 
of practical importance, it has been noted that in situ sediment 
shearing resistance versus rotation data at different rotation 
rates can be obtained using the wire line vane shear device of 
McClelland Engineers. It was thus deemed desirable to develop 
a laboratory vane shear test and analysis procedure which could 
be evaluated and, if valid, extended for use with in situ data. 
This section derives the equations and outlines the computational] 
method used. The theory presented here is an exposition of the 
proposed procedure of Schapery (34) for computation of the visco- 
elastic shear modulus from vane test data. The procedure involves 
initial solution under the assumption of linearly elastic sedi- 
ment, from which the linear viscoelasticity solution can be 
obtained; and finally a graphical characterization of the nonlin- 


earity of the sediment is derived. 





be: 


Calculation of Shear Modulus from Vane Test Data - Elastic 
Assumption. - First consider the sediment to be an infinite mass 
of homogeneous, isotropic linearly elastic clay. Consider an 
embedded vane rotated through an angle 9. Assuming that the vane 
is rotated fast enough to insure undrained conditions and migra- 
tion of pore water is not important, the applied torque, T, may 
be characterized by: 

Hiemttch We C10) coe Gea a Wl a Gal (6) 
where a and L are the vane radius and length, respectively, and 
G is the elastic shear modulus. Eq. 6 can be alternately expressed 
as: 

OF- Salsa, Las G, 0) 
which, by dimensional analyses and the Buckingham Pi Theorem, 


can be expressed as: 


and can be solved for T: 
ey rs ek 

Te Gaol f(s 6) 
The linearity assumption implies that T is proportional to 6. 

Applying this to Eq. 9 and rearranging the geometry terms yields: 

_ a 

T = Gea°L-f, (7) 


The quantity ae 1s a constant which depends at most on dimen- 


Sionless ratios defining the vane geometry. It is noted that this 





result was obtained using only the vane length and width as defin- 
ing the vane geometry. If other geometry factors are used, how- 
ever, the same result is obtained, i.e., f, is a dimensionless 
constant for any vane geometry, and depends on geometry ratios 
only. If f, is defined as the vane geometry factor, B, then 
Eq. 10 becomes: 
T = Ba*LGe 
Since B depends only on geometrical ratios it will be the same 
for all geometrically similar vanes, regardless of absolute size. 
(This fact will be used later to extend the laboratory study to 
the in situ case.) 
Under the elastic assumption Eq. 11 can be used to compute 


the modulus G if torque versus rotation data were obtained. It 


is first necessary, however, to determine the numerical value of 
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the vane geometry factor B. Methods of determining B are described 


at the end of this section. Attention will now be turned to the 
viscoelastic case. 

Obtaining the Linear Wiscoclasticit? Solution from the 
Elastic Equation. - The above discussion has been limited to 
characterization of a linearly elastic material. £&q. 11 consti- 
tutes the elastic solution for a rotating vane in the clay mass. 
For many viscoelasticity problems and the present one in partic- 
ular, the linear viscoelastic solution can be found from the 


Laplace Transform of the elastic equation, so long as the kind 
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of data (strain, strain rate, etc.) specified at any specific 
boundary point is the same at all times. This statement is known 
as the correspondence principle (30). The correspondence princi- 
ole is valid as long as the material response can be described 

by constitutive equations of the convolution form which implies 

the material is nonaging and temperature dependency of the response 
is not important. Each of these conditions being satisified in 

the case of the vane shearing of clay, the correspondence prin- 
ciple may be applied. First, taking the Laplace Transform, 
L{f(t)}, of each of the variables in Eq. 11 yields: 

CEO Loa Olt) Gauss. : 20k. ee we (12) 
Then, applying the correspondence principle the linear viscoelastic 
solution can be obtained by replacing G with s 2{G(t)}, where 
G(t) is the relaxation modulus in shear, and s is the Laplace 


Transform parameter, and finding the inverse. 


le Babee. 1G} 7{0) aeemeeoih, st io 6) ci Me ens (3) 
Taking the inverse transform yields: 
C de 
Ts Batt | Gtr) Pdr... ee ee, (14) 


0 

where G(t) is the viscoelastic relaxation modulus in shear, t is 

a dummy time variable of integration and (t-t) represents the time 
elapsed since rotation began. Eq. 14, then, is the viscoelastic 
equivalent of Eq. 11 and we see that the result amounts to the 


replacement of the product GO with the so called convolution 





integral, which forms the basis for linear viscoelasticity theory. 
For constant rotation rates Eq. 14 can be simplified into a 


form analogous to Eq. 11 by first letting 


where C is a constant and thus, at any specified time t: 


- 
eee ees ee ere (16) 


Then d6/dt can be removed from the integral to obtain: 


t 
T= Batt & | asa 0) ool eee (17) 
0 
dé _ 
for RT C 


Further, a change of variable in the integral: 


t' = (t-t) at T#= 0, t'et 
dt' = -<dt Gah ae 
yields: 
9 t 
T = Batt & | Git’) dee... oo. =. (18) 
0 
or 
1 t 
T = Ba’ @ [> | att’) dt! a as a (19) 
) 


It is recognized that the term in brackets in Eq. 19 is the time 
averaged value of the shear modulus G(t) and denoting this by G, 
Eq. 20 is obtained: 

fmeeGa leno Glo. eee 2 ek ee ee te ke (20) 
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(es, 


where 
_ 1 t 
G = G(t) = | GUCCI Se es meet 2) 
0 


Now the viscoelastic relaxation modulus of clays can often 
be characterized as a power Jaw in time (25, 33). If such a rela- 


tion is assumed for G(t) we obtain: 


where G, and n are constants. Using the definition of G above 


(Eq. 21) we obtain 


G= | G, 417" GU MCR. weCw eet, 6 et. ee (23) 
0 
or 
Ga : 
a -N 
G i T-n t « ee e® e© e8® e© e@ e# #* &® e® e@ 8 &© © 8 @ @ «@ ‘.* @ ° (24) 


Evaluation of G, and n in Eq. 22 would constitute a solution to 
the linear viscoelasticity problem. Eqs. 20 and 24, together, 
provide the framework for determining these constants from con- 
stant rate vane shear data. The computational! procedure is out- 
lined below. 

First, however, it is of interest to note that in applications 
where sediment response to cyclic or vibrational inputs is desired 
(such as the wave induced slope stability problem) then a fre- 
quency dependent complex shear modulus, G*, is required, where: 


GSE GUE ie lckvee en Seeeeets en 5 2 SRG cs SM RES. Goi. 1s 0 et (25) 
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G' and G" are constants and i = v-1. G* can easily be obtained, 
given the input frequency, w, and the values of G, and n of Eq. 


22. G' and G" are computed from: 


G = [GT(i-n) cose] wo 2... ee (26) 
a= (G,r(len) sin da ww ww (27) 


where T(1-n) is the Gamma function with argument I-n. Eqs. 26 and 
27 are obtained from well known viscoelastic theory, given the 
power function, Eq. 22, and a sinusoidal input. 

Computational Procedure and Characterization of Non- 
linearity. - It is now shown how Eqs. 20 and 24 can be used to 
characterize the material from vane test data. Suppose a series 
of vane tests were run on the same material, each test at a dif- 
ferent constant rotation rate, 0 = &. Continuous measurement 


of torque versus rotation angle would yield data such as shown 


in Fig. 8, where the ordinate variable is normalized torque, T,: 


By selecting an angle, say 6,, G can be computed for each rotation 
rate Q.> Oo» 0, etc. by Eq. 20, j.e.: 
Eco eyes" all eer fats ule Bielktet alle > use) -« «© sl "s « “se « (29) 
Ba“L6@ 98 


By determining the corresponding time for each 6, the plot of log 


G versus log t shown in Fig. 9 for @, can be obtained. 
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FIG. 8. - NORMALIZED TORQUE VERSUS ROTATION ANGLE 
AT VARIOUS ROTATION RATES 
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FIG. 9. - MODULUS VERSUS TIME FOR VARIOUS 
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Now, if the real behavior of the material can be represented 
by the power law (Eq. 22) then from Eq. 24 we see that the log G 
versus lot t plot at 6, will be a straight line of slope -n. 

Also, if the material is linearly viscoelastic, by definition 

(see Eq. 21) G is not dependent upon 6 and the log G versus log 

t plot for all angles will be the same straight line. Significant 
departure from the assumed power law can be noted if the deta do 
not plot on a straight line, and nonlinearity (6 dependence) is 
revealed by separation of the plots for different angles 6,> 5,» 
6. etc. as shown in Fig. 9. Thus direct assessment of the non- 
linearity is available from examination of the graph. 

Once Fig. 9 is obtained n is found from the slope (slope = -n) 
and G, is obtained by selecting any time, finding the correspond- 
ing G from the figure (at the desired 6) and computing G, from 
Eq. 24, 1i.e.: 


With n and G, determined G(t) can be determined for any t using 
Eq. 22, and the complex modulus, G*, can be determined using Eas. 
25 through 27. 

Determination of the Vane Geometry Factor, B. - In the dis- 
cussion above the factor B was obtained as a dimensiontess con- 
stant for any particular vane geometry, regardless of absolute 


size. In order to use vane shear data to find G(t) it is 





uit) 
=> 
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necessary to evaluate B for the vane shape being used. 6B can be 
determined in several ways: 

1. By theoretical estimate. 

2. By running vane tests on a material of known viscoelastic 
properties and computing B from Eq. 20. 

3. By running vane tests on a thick fluid of known viscosity 
and computing B from the relation T = Ba2Lon where n is 
the viscosity of the fluid. 

4. By measuring G for the sediment by some other means and 
using it, with vane test data, in Eq. 20. 

These methods will be explored further here. 

It is first useful to derive the theoretical estimate. To 
accomplish this consider a homogeneous, isotropic, clay mass of 
infinite extent. Replace the vane by a long rigid cylinder of 
radius, a. The clay is considered linearly elastic and is assumed 
to adhere to the surface of the cylinder. As the cylinder is 
rotated the shear stress across the cylindrical surface is uni- 
form. The shearing stresses and deformations in the clay sur- 
rounding the cylinder can now be derived in the manner presented 
by Cadling and Odenstad (4). 

Fig. 10 depicts the shearing stresses acting on a volume 
element a distance, r, from the center of rotation. Due to the 


symmetry of rotation the element is subjected only to the net 


— 





34 


shearing stresses t, and on dt. The equation of equilibrium 


for a radial slice is: 


ca a Goo 6 Sao 6 co 5 i oer (31) 
where t, is the shearing stress at the cylinder surface. Now 
introduce the relation: 

i ee So ME ow Ss (32) 
where G, as usual, is the elastic shear modulus and if is the 
shearing enaati at r. Combining Eqs. 31 and 32 yields: 

Vp 222 Yq I awe... (33) 


where y, is the shearing strain at the cylinder surface. As 
depicted in Fig. 11, during rotation each clay particle moves a 
short distance on a circle of radius, r. The particles at any 


r move through an angle: 


If the particles at r + dr were to move a distance (r + dr) 6 on 
a circle, then no shearing would occur in the element at r. In 
reality, the particles at r + dr move a distance (r + dr)(6 + dé). 


Hence, the strain at r is: 


~ (Veemcles) (alae dr (8 + do) (35) 
ie ae se TO igie Sle) 6 


Yoo -(r + dr) Ais age ©) 0 0 UNO c so ncn ce (36) 


which, as the size of the element is allowed to vanish, yields: 


dé 
iweeOady me ts ee ee ee (37) 








FIG. 10. - SHEARING STRESSES ACTING ON ELEMENT 
OF CLAY DUE TO ROTATING CYLINDER 





FIG. 11. - DEFORMATION OF CLAY SURROUNDING 
ROTATING CYLINDER 
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Now, combining this result with Eq. 33 and dropping the negative 


sign we obtain: 


where C is determined when 6 = 0 at r = ~ and is, C = OQ. There- 


fore: 


With this-result in mind we now turn to the determination of the 
torque on a rectangular vane of length L, in terms of the shear- 
ing stress t, at the surface of the cylinder of revolution (the 
assumed failure surface). Neglecting end effects, the torque is 
the moment. produced by the shearing stress T, on the surface of 
the cylinder, which is: 

UC CCl Ee Cee ee nee (43) 
substituting t, from Eq. 42 yields: 

Game a Lede, G 
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and since 0, is just the angle through which the vane is turned 
we now drop the subscript and compare Eq. 44 with the Eq. 11 to 


obtain: 


It is emphasized that this is merely an estimate. In practice B 
must be determined from one of the empirical procedures listed 
above, especially if the vane blades are not rectangular and the 
surface of revolution of the vane is not a cylinder. 

Pursuing this development further, however, will yield an 
estimate of the effect of not having an infinitely wide mass of 
clay, i.e., the effect of container size. Assume that the vane 
is rotated about the longitudinal axis of a cylindrical sediment 
sample in a container of inside radius, b. The shearing displace- 
ment at the inside of the container will be zero. Making use of 
Eq. 39 but this time evaluating the constant, C, with @ = 0 at 
r = b yields: 


C 


] 1 
Q@ = oF a* (or ~ BT) cS a 5. ee . (47) 


(sae BO (48) 
an (pe > pe) 
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Combining this with Eq. 48 at r = a, solving for T and rearrang- 


ing, yields: 


2 
Oa (50) 
(1 - #2) 
Comparing this with Eq. 1] yields 
B = at See, ee, ek a (51) 
(1 - Bz) 


From this result we see that if a is less than about .3b then the 
ratio a7/b* will be small, B = 4m, and the effect of the container 
will be negligible for al] practical purposes. 

Attention is now turned to the empirical methods for deter- 
mining B. The most convenient, from an experimental standpoint, 
are methods 2 and 3 (see page 33). However, method 2 is imprac- 
tical because it is difficult to find linearly viscoelastic 
materials which can be used as standards. Method 3, based on the 
equation: 

T = Ba*Lné 
is more easily applicable, where n is the viscosity of a reference 
fluid. This equation is obtained in the same manner as was Eq. 

11, where the linearily elastic material is replaced by a Newtonian 
fluid of viscosity n. For experimental purposes the fluid used 
should offer shearing resistances the same order of magnitude as 


those expected from the clay sediments at the rotation rates to 
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be used. The use of this technique is demonstrated in the follow- 
ing section on experimental procedure. 

The last method available for determining B is to measure G 
of the sediment by an independent experimental method and compute 
B using G, vane test data and Eq. 20. One such independent method 
for finding G is measurement of sediment response to a vibrational 
input. Briar and Bills (3) have demonstrated such measurements 
in solid rocket propellants. Their method employs a small peizo- 
electric disc which is inserted in the sample and electrically 
excited. The voltage and frequency response output, together 
with the response obtained in air can be used to compute the 
complex modulus, from which G can be easily obtained (3). The 
technique has also been used by Cohen (6) and others, on submarine 
sediment samples. Cohen's device was a piezoelectric bar excited 


into torsional vibration. 
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EXPERIMENTAL PROGRAM 


The laboratory phase of this research involved development of 
a precision laboratory vane shear device and test of submarine core 
samples at different rotation rates. As discussed below, the 
samples were obtained from three boreholes off the coast of 
Louisiana. The wire line vane shear device developed by McCleiland 
Engineers was used to measure torque versus rotation, at several 
rotation rates, and at several locations in one of the boreholes. 
On the basis of the theoretical development above, the laboratory 
vane was constructed to be geometrically similar to the wire line 
vane, thus making possible a comparison of the laboratory and in 
situ rate tests. 

The laboratory vane shear apparatus was designed to meet the 

following primary criteria. 

1. Continuously record torque versus rotation data while 
rotating a vane at constant rotation rates in sediment 
core samples. 

2. Operate at rotation rates ranging over several decades, 
with the rates available in the McClelland wire line 
vane included. 

3. Incorporate a vane geometrically similar to that used 
in the McClelland wire line device to make possible com- 


parison of results. 
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Description and Operation of Apparatus. - Fig. 12 shows the 


experimental apparatus for the vane tests, and Fig. 13 shows the 
details of the vane shear device. The device consists of a small 
variable speed D. C. motor which drives the vane shaft through a 
Step function speed reducer with output to input speed ratios of 
1:1 to 1:1000 in 10 steps. The total rotation speed range thus 
obtained is 5x10"° radians/sec to 0.57 radians/sec. Torque on the 
vane due to sample shearing resistance is measured by a strain 
gaged torque transducer and is recorded on a Visicorder recording 
oscillograph as a plot of torque versus time. The torque sensi- 
tivity of the instrument is 0.001 in.-Ilb, with full scale settings 
from 0.2 in.-1b to 10 in.-1b available. With the vane and motor 
used, stiffness data on sediments ranging in stiffness from less 
than 10 psf to over 500 psf can be obtained. Rotation is measured 
by the change in resistance of a single turn linear, precision 
potentiometer connected to the vane shaft by anti-backlash gears. 
The gears provide a vane shaft to potentiometer shaft rotation 
ratio of 1:3.67. Rotation readout is a continuous plot of rota- 
tion angle versus time on a second channel of the Visicorder. 

Vane rotation changes of less than 0.005 radians are resolvable. 
Vane operation is controlled via a control module which houses 

the on-off-reverse switches for the vane motor and circuitry 
associated with the rotation transducer. The vane is shown in 


Fig. 14. 








FIG. 12. - TEST APPARATUS FOR LABORATORY VANE 
SHEAR RATE TESTS 
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FIG. 


13. = DETAILS OF VANE SHEAR DEVICE 
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FIG. 14. - THE TEST VANE 
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The apparatus includes a pulley arrangement and torque trans- 
ducer calibration jig which facilitates routine physical check 
calibrations of the torque transducer. By simply disconnecting 
a shaft coupling between the upper vane shaft and the step func- 
tion speed reducer, and disengaging the gears to the rotation 
transducer, the upper end of the shaft is free to turn in a pair 
of precision ball bearings. By restraining the lower end of the 
transducer in the calibration jig and applying torque to the upper 
end via the pulleys by a string and weights, calibration is 
obtained. The calibration configuration is shown in Fig. 15. 

The vane shear device is mounted on the uprights of a Wykeham 
Farrance Soil Test Machine. The machine's base pedestal serves 
as the sample pedestal and provides a means to very smoothly 
push the vane into the sample by slowly raising the pedestal. 
Detailed descriptions of the apparatus components and a descrip- 
tion of initial calibration and check out procedures are found 
in Appendixes III and IV. 

Sample Selection. - As shown in Fig. 12 (page 42) the test 
samples were 5 in. long, 2.125 in. diameter core tube samples. 
They were recovered from three boreholes (B-1, B-2, B-3) in the 
south Pass area of the Mississippi Delta off the Louisiana 
coast. The core samples were obtained by the down hole wire line 
core sampler and in situ vane shear rate tests were conducted at 


several locations. The in situ vane tests were accomplished by 
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FIG. 15. - TORQUE TRANSDUCER CALIBRATION ARRANGEMENT 
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McClelland Engineers, Inc. using the wire line vane apparatus 
previously discussed (9). After recovery the samples were sealed 
and transported to the Texas A&M University campus, where they 
were stored in a controlled humidity cool room until used. 
Samples from the sediment depths shown in Table |] were 


selected for testing. These samples were selected to yield 


TABLE 1. = SAMPLES TESTED 


ee ee ee 


Borehole B-1 Borehole B-2 Borehole B-3 


Ce Oma 55 ft 
40.5 ft 


(4B) DE 
104.5 ft 





representative sediment profiles for boreholes B-1 and B-2 and 

to sample an area of extremely low strength in B-3. In situ rate 
test data were also available near each of the depths selected in 
borehole B-1. 

An additional sample selection criterion was used. This was 
the degree of sampling disturbance as evidenced by X-rays. Fig. 
16 is an X-ray photograph of an acceptable sample. Fig. 17 shows 
a sample that was rejected on the basis of excessive fracturing 
and nonhomogeneity, probably caused by gas expansion on recovery. 


Appendix V presents X-rays of all samples tested in this study. 
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FIG. 16. = X-RAY PHOTOGRAPH OF ACCEPTABLE SAMPLE 





FIG. 17. = X-RAY PHOTOGRAPH OF UNACCEPTABLE 
SAMPLE SHOWING FISSURES DUE TO GAS EXPANSION 
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Experimental Determination of Vane Geometry Factor - B. - The 


vane geometry factor for the test vane configuration was determined 
by rotating the vane in a reference viscosity fluid and computing 
B from Eq. 52. The fluid was Dimethylpolysiloxane silicon fluid 
sold under the trade name Dow Corning 200. Its rated viscosity 
at 25°C is 5000 poise. The fluid was calibrated at Texas A&M 
in a Weissenberg Model R-17 Rheogoniometer and found to be 
Newtonian over the required shear rate range with a viscosity of 
5507 poise at 20°C. The vane tests were conducted at 20°C. Turn- 
ing the vane in this fluid resulted in torques at high rotation 
rates roughly equivalent to those encountered in very soft sub- 
marine sediments. 

The value of B determined by this method was 9.596. It is 
noted that this is somewhat below the theoretical value of 
4n = 12.56. This large difference is due to the fact that the 
vane does not have rectangular blades (see Fig. 14, page 44), 
as assumed when the theoretical value was determined. Thus the 
use of the overall vane length L, in Eq. 52, results in a smaller 
value of B. 

sample Test Procedure. - Sample testing was accomplished by 
running vane tests at several rotation rates at adjacent locations 
“mali Sample. The four rates were, nominally 0.26, 0.026, 0.013, 
and 0.0013 radians/sec. The two middle rates were selected to 


approximate the rates used in the in situ wire line device, and 
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the high and low rates provided the wide range deemed necessary 
to obtain sufficient data for use of the computational procedure 
discussed above. Continuous recordings of torque and rotation 
angle versus time were obtained for each rotation rate. During 
the testing sequence for each sample the samples were resealed 
and replaced in the cool, humid room whenever they were not 
actually being used. This eliminated moisture content changes. 
The test sequence on each sample was as follows. 

1. Upon being opened each sample was visually inspected 
and a conventional vane test of ultimate strength was 
run using a motorized Wykeham Farrance miniature vane 
(Fig. 18). This test was run to obtain index lab 
strength values for the samples and to estimate the 
range of torques expected in the rate tests so the 
appropriate full scale torque range could be set on 
the recorder. 

2. A physical calibration of the torque transducer was 
conducted over the range of expected torques. 

3. Appropriate motor voltage and step function speed 
reducer settings were made to obtain the desired rota- 
tion rate, and the vane was rotated to obtain a zero 
rotation reading on the recorder. 

4. A test run was then began by holding the sample on the 


pedestal, raising the pedestal, and forcing the vane 
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FIG. 18. - WYKEHAM FARRANCE MOTORIZED MINIATURE 
LABORATORY VANE SHEAR DEVICE 





into the sample at an undisturbed location. The recorder 
drive and vane motor were then actuated and the torque- 
rotation-time data were obtained at that rotation rate. 
Fig. 19 is an example of data so obtained. 

5. Upon completion of each test, at one rate, the sample 
was removed and Steps 3 and 4 were repeated at the next 
rotation rate. 

6. Immediately after the vane tests were complete moisture 
contents were obtained for each sample. In addition 
Atterberg limits were later obtained for each sample. 

Data Reduction. - Each vane test on each sample resulted in 

data such as shown in Fig. 19. The data were reduced by selecting 
an angle, using the rotation calibration (see Appendix IV) and the 
known time scale on the chart to determine the corresponding time, 
and, using the torque transducer calibration data obtained for 
each set of runs, the corresponding torque was measured. Such 
torque-rotation-time data were obtained at 0.0] radian increments 
from 0 to 0.1 radians and at 0.15 and 0.20 radians for each 

run. Maximum torques (failure) normally occurred at or below 

0.20 radians. The rotation rate for each test was not determined 
from instrument settings but was computed as the average rate to 
an angle of 0.10 radians computed from the angle-time data. 

Since the actual time associated with each angle, not the rotation 


rate, 1S used in the modulus calculations, this approximation of 


De 
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rotation rate is of no consequence. The measured rotation angles 
were corrected for shaft twist between the vane and the rotation 
sensor using a correction factor calculated as the elastic shaft 
twist at the measured torque. The corrections were slight. 
Torque, T, versus rotation angle, 6, plots were then made 
from the raw data for each sample, as shown in Fig. 20. Appendix 
VI contains these plots for all samples. From these plots the 
normalized torque, T,, versus 6 data were obtained, the values 
of the time averaged viscoelastic shear modulus, G, were found, 
and G, and n of the power law (Eq. 22) were determined using the 
log-log plotting method shown in Fig. 9, page 31. Fig. 21 shows 
the log-log plot for Sample B-1 at 46 ft. Plots for all samples 


are in Appendix VI. 
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RESULTS AND DISCUSSION 


Application of viscoelastic theory and the power law, Eq. 22, 
to vane shear data have yielded values of n and G, for each of the 
test samples. Examination of the log-log plots of Appendix VI 
indicates that n is very nearly independent of 8. However, due 
to the nonlinearity of the sediment, 

Gh iaG (0) OM SEES cot? Sipe. 3 om (52) 
and thus the log-log plots at different angles do not fall on the 
same line. Values of G, were computed at each of the several 
angles used for constructing the log-log plots. In addition, an 
estimate of an initial G,, G, at 6 = 0, was obtained by plotting 
G, versus 6 for each sample on linear paper and extrapolating the 
graph to 6 = 0 using a straight line from the lowest 6 data point. 
This procedure is common practice in viscoelastic theory to esti- 
mate the modulus at very short times from experimental data. 

Experimental Results. - Table 2 presents the experimentally 
obtained values of n and G, for each sample. 

The plots of G, versus 8, showing the angle dependency of G,, 
are presented in Figs. 22 through 31. Review of Table 2 and these 
figures indicate, as would be expected, a general decrease inn 
with depth and an increase in G,. This indicates decreasing time 


dependence of the modulus with depth and an accompanying increase 
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in stiffness. Fig. 32, a plot of n and G, at 6 = O versus depth, 
graphically depicts these variations with depth. It is noted that 
for sample B-3 at 55 ft the n and &, data appear anomalous. How- 
ever, as was previously noted, this sample was from an extremely 
weak region. 

The conventional laboratory soil mechanics tests on the 
samples are summarized in Table 3. It may prove useful to con- 
sider the possibility of using conventional parameters, such as 
liquidity index, L.I., and ultimate shear strength, Cr as pre- 
dictors of G, and n for use with power law modulus computations. 
Figs. 33 through 36 are plots of G, and n versus L.I. and C,. 

The available data indicate strong relationships between L.I. 
and n, and between Cy and G,. This is consistent with the con- 
cept of viscoelasticity in clays, especially the correlation of 
n, the indicator of time dependence in the power law, and L.I., 
an indicator of the water content relative to the water content 
at which the clay tends to act as a fluid. At this time the 
available data are insufficient to attempt quantification of 
these relationships. 

Comparison of Laboratory and In Situ Results. - Limited in 
situ torque versus rotation angle data, at several locations jn 
Borehole B-1, were obtained during the original sampling. Data 
were obtained at three rotation rates. Three of these tests were 


selected for comparison with the laboratory results since they 
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were conducted at depths corresponding with laboratory sample 
depths (46 ft, 61 ft and 100 ft). 

It is noted that the in situ tests were conducted at rota- 
tion rates between 0.0182 rad/sec and 0.0077 rad/sec. This 
narrow range (as opposed to a range of about 10° rad/sec for the 
laboratory tests) precludes the use of these data to compute values 
of G, and n as was done with the laboratory data. Time differ- 
ences in the points obtained for the log G - log t plots are so 
small that the slope of the plot changes greatly for small changes 
in G. Thus experimental error in torque, angle and time measure- 
ments would be amplified. For these reasons the in situ rate 
test data are not considered sufficiently reliable to attempt 
direct comparison of laboratory and in situ values of G, and n. 

It is possible, however, to obtain a check on the validity 
of the power law representation of the shear modulus by using 
laboratory values of &, and n to predict the in situ torque versus 
rotation data. This can be done by considering the normalized 


torque: 


which can easily be computed for each in situ data point. Due to 
the vane similarity this can be directly compared to T, predicted 
from the laboratory results. To do this substitute the power law 


for G(t) into Eq. 14 to obtain: 
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me - 
"> Batt * J G, (t-1 


Now, it is noted that G, is dependent on 8 but is not time depen- 
dent. If the rotation rate is constant then 2 C. By using 


this and carrying out the integration: 


G (9) 1l-n 
Lied ] 
i, Se 


Now since @ = C then at any time, t, 6 = Ct and this can be sub- 


Stituted into Eq. 54 to yield: 


: G, (9) att 
-n 


Th 
Now, if it is desired to predict the in situ T, at a particular 
6 - t pair (i.e. at angle, 6, with rotation rate 2) Eq. 55 can 
be utilized by using n determined from a laboratory sample from 
the location where the prediction is desired and using the labora- 
tory value of G, corresponding to the 6 of interest. 

These predictions have been calculated to angles of .06 radi- 
ans for the Borehole B-1, 46 ft, 61 ft and 100 ft depths. Figs. 
37, 38 and 39 show the comparisons. The predictions could have 
been carried out to larger angles if desired. In light of the 
fact that the in situ shearing resistances are expected to be 
Slightly higher than those of samples (28), the comparisons are 
remarkably close. The large spread of some of the in situ data 
from the predicted curves is believed to be primarily due to non- 


homogeneity of the sediment and certain uncertainties associated 
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with the in situ tests. The comparisons tend to verify the proce- 
dure for determining the vane geometry factor and lend experimental 
evidence to the theoretical assertion that it is independent of 
absolute vane size, and that pore water migration is not important 
in this clay at the vane rate used. They also support practical 
use of the power law representation of the time dependent shear 
modulus of submarine sediments. 

Comment. - The experimental results appear encouraging. The 
power law viscoelastic shear modulus and the vane test and analy- 
Sis technique seem to lend themselves to practical determination 
of the modulus and prediction of in situ sediment stress-deforma- 
tion behavior. The simple form of the equation makes it easily 
applicable in many types of problems including response to vibra- 
tional inputs. Available data indicate the procedure yields 
accurate results. However, to verify the validity, additional 
testing is necessary. The general procedure (ie. application of 
the power law) might be verified by determining G(t) by some other 
means and comparing with G(t) predicted from the vane data and the 
power law. One other means might be the use of the vibrating disc 
technique (3). The equality of the vane factor B for the in situ 
and Jab vanes might be further verified by running the two vanes 
in the same homogeneous clay material in the laboratory environment. 

It is noted that no direct measurement of the in situ shear 


modulus was deemed possible in this study. This is only a result 





he 


of inadequate rotation rate range and uncertainties in the in situ 
angle-torque-time determinations. Development of an in situ vane 
device which can rotate at speeds about three decades apart and 
employ more accurate means to determine actual time versus rota- 


tion and torque, continuously, is required. 
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CONCLUSIONS 


Submarine clays exhibit viscoelastic stress-deformation 


characteristics. This study has presented a practical method 


for determining a viScoelastic shear modulus for the sediments. 


On the basis of the theoretical and experimental considerations 


it is concluded that: 


I 


The time dependent shear modulus of submarine sedi- 
ments can be reasonably represented by the power law 
equation, G(t) = GRetier, where G, is strain dependent. 
The vane shear test, run at several rotation rates 

and measuring continuous torque-rotation-time data, 
provides a practical means of determining the shear 
modulus, both in the laboratory and in situ. 
Laboratory vane shear data can be used to estimate 

the in situ stress-strain behavior, and, if sampling 
disturbance is neglected, to predict the in situ 
modulus. 

The presently available wire line down hole vane shear 
device of McClelland Engineers does not have a suffi- 
ciently large range of rotation rates to utilize it 
for direct in situ modulus determinations. Also, 
refinement of the measurement methods for torque-rotation- 


time data is required for use of this instrument. 





ay 


6. 


The experimental data indicate a relationship between 
liquidity index and the constant n in the power law. 
The experimental data indicate a relationship between 


the maximum sediment shear strength and the initial 


(6 = 0) value of G,. 
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APPENDIX II. - NOTATION 


vane radius 

vane geometry factor 

a constant 

maximum vane shear strength 

vane diameter 

viscosity 

elastic shear modulus 

viscoelastic relaxation modulus in shear 
average viscoelastic modulus in shear 

a constant in time 

Shear strain 

Shear strain rate 

vane length 

liquidity index 

Laplace transform of f(t) 

a constant 

frequency 

plasticity index 

radial distance from center of rotation 
torque 


normalized torque 
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Mt] 


Wn 


time 

dummy time integration variable 
shear stress 

rate of shear stress 


angle of vane rotation 


= liquid limit 


= natural water content 


plastic limit 
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APPENDIX III. - DESCRIPTION OF EQUIPMENT 


This Appendix is intended to provide a detailed description 


of the equipment. The test apparatus and a close up of the vane 


shear device are shown inFigs. 12 and 13 (pages 42 and 43). 


The components and arrangement are listed and described in detail 


below. 


List of Components 


1. Vane Shear Device 


(g) 
(h) 


D.C. Motor 

Step Function Speed Reducer 
Bearing Plate and Gear Assembly 
Torque Transducer 

Vane Rotation Transducer 

Vane 

Mounting Assembly 


Calibration Assembly 


2. Stand and Sample Pedestal 


3. Motor Control and Transducer Circuity Unit 


4. Torque Transducer Signla Conditioner 


5. Motor Power Supply 


§. Strip Chart Recorder 
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Description of Components 


ie 


Vane shear device (see Fig. 13, page 43). 

(a) D.C. Motor - The vane rotation motor is an American 
Electronics Inc., No. 3255P reversible D.C. gear motor. 
It is operable between 5 and 35 VDC and, incorporating 

a gear reduction of 4126:1, produces output rotation 
ORES of 0.5 to 5.5 rpm. Continuous duty torque is i0 
in-oz. Motor dimensions are 1.375 in. diameter by 3.31 
in. long. The motor is mounted to the back mounting 
plate with an aluminum bracket and vibration isolating 
rubber grommets. The output shaft is connected to the 
input shaft of the speed reducer with 0.25 in. stainless 
steel shafting and a vibration reducing semi-flexible 
coupling. 

(b) Step function speed reducer - The speed reducer is 
an INSCQ Model No. 00740 Step Function Speed Reducer 
which provides rotation output to input ratios of 1:1, 
ecg le Cee lO, bacO, 1:50, 03:t00, 1:200, 17500 and 
1:1000. Changing the ratios is easily accomplished by 
turning a knob on top of the unit. The input and output 
shafts both come out of the bottom of the unit and are 
0.25 in steel. The reducer is rated at 100 in-oz torque. 
The reducer is attached directly to the back mounting 


plate by machine screws. 
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(c) Bearing Plate and Gear Assembly - The bearing plate 


and gear assembly is mounted to a mounting plate on the 
bottom of the speed reducer by four machine screws with 
0.9 in. sleeve stand offs. The Assembly consists of 

two 3/10 in. thick steel plates which are held together 
by four 2 in. spacer studs. The plates provide mountings 
for four precision, 0.25 in I.D., miniature ball bearings; 
two for the upper end of the vane shaft (between the 
torque transducer and the speed reducer) and two for the 
rotation transducer shaft. Qn the shafts, between the 
two bearing plates, are mounted the gears for the rota- 
tion measurement and the 1.5 in diameter torque trans- 
ducer calibration pulley. The gears are anti-back lash 
stainless steel spur gears providing a rotation step 

up from vane shaft to rotation transducer shaft of 
1:3.67. The bearing plates are slotted on the motor 
side to accomodate the motor to speed reducer shaft. 

(d) Torque Transducer - The torque transducer is mounted 
as part of the vane shaft just below the bearing plate 
assembly. It is 2024-T3 aluminum 0.75 inch bar stock, 
2.25 inches long with a 0.250 in. diameter center drilled 
hole running its length. The middle portion of the 
transducer is necked to provide 0.015 inch wall thick- 


ness and design stresses of 6000 psi in the outer fiber 





at 150 in-oz torque. Torsional strain in the transducer 
is detected by a four arm balanced strain gage bridge. 
The strain gages are SR-4, 350 ohm. The bridge is con- 
nected directly to the signal conditioner. 

(e) Vane Rotation Transducer - The vane rotation trans- 
ducer is a Beckman Helipot Model 6103, 10,000 ohm ball 
bearing precision potentiometer. Full scale linearity 
is rated at +0.5%. The potentiometer is mounted to the 
back mounting plate with an aluminum bracket. The 
potentiometer shaft is driven by the vane shaft through 
the two anti back lash gears in the bearing plate and 
gear assembly. A 5 volt reference potential from the 


control unit is across the outside terminals of the 


potentiometer. Rotation is measured by measuring voltage 


changes at the rotor terminal as the potentiometer (and 
vane) rotates. 

(f) Vane - The test vane is Shown in Fig. 14, (page 44) 
and is composed of four stainless steel blades silver 
soldered onto a 7 in. length of 0.25 inch stainless 
shafting necked to 0.125 in. at the vane end. The vane 
blades are 0.0156 in. thick. The total length of the 
vane blades (tip to tip), L, is 1.15 in. and the total 
width is 0.50 in. (i.e.-a = 0.25 in.). Both ends of 


the blades are tappered at 45° angles. This vane is 
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geometrically similar to the vane used in the McClelland 
Engineers wire iine vane shear deivce. The vane is 
connected to the bottom of the torque transducer by a 
shaft coupling. 

(g) Mounting Assembly - The mounting assembly provides 
the structure for mounting the vane shear device to the 
uprights of the Wykeham Farrance Soil Test Machine. The 
assembly includes two horizontally mounted 0.5 in. steel 
top and bottom plates with holes in each end such that 


they slide over the uprights. (see Fig. 12, page 42). 


These plates are held in place an the uprights by collars. 


Between the top and bottom plates the .375 in. thick 
steel back plate is mounted vertically. The other com- 
ponents of the vane shear device are mounted to this 
plate. 

(h) Calibration Assembly - The torque transducer 
calibration assembly is illustrated in Fig. 15, page 
46). It is composed of the 1.5 in pulley mounted on 
the vane shaft, a second pulley mounted as shown in 
Figure 15, a piece of fishing line, calibration weights, 
and a three point clamping jig for immobilizing the 
bottom end of the torque transducer. The assembly is 
used by first removing the vane and disconnecting, and 


Sliding up, the shaft coupling between the upper vane 
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shaft and the step function speed reducer. The fishing 
line is then connected to the vane shaft pulley and led 
Over the other pulley, and the lower end of the torque 
transducer is retained from rotation by the clamping 

jig. The restraint is applied by a carefully tightened 
three set screw jig to avoided putting bending moments 

on the torque transducer. Mercury filled weights cali- 
brated in even in-oz increments (computed on basis of 1.5 
in. diameter pulley on vane shaft) are then hung from the 
line and the torque readout is obtained on the Strip 

chart recorder. 

Stand and Sample Pedestal - The uprights and pedestal of a 
Wykeham Farrance Soil Test Machine, Model T-56-B are used 
to support the vane shear device and provides a means to 
raise the sample, pushing it onto the vane. ( SCCRiaiGn 

12, page 42). The pedestal is operated by a hand crank. 
Motor Control and Transducer Circuity Unit - This unit 

was assembled to provide a convenient means to Start, stop 
and reverse the vane motor and to Provide a constant volt- 
age source and calibration system for the rotation 
transducer. Fig. 40 shows the unit and Fig. 41 is a 
schematic. The motor control circuit consists of an on- 
off push button weitch and a motor voltage polarity 


reversing switch. The rotation transducer circuit is 








FIG, 40. - MOTOR CONTROL AND TRANSDUCER CIRCUITRY 
UNIT 
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driven by 28 VDC from an external battery. This is 
reduced to 5 volts by the reference voitage transformer. 
When the control switch is in the "run" position the 5 
volts are applied across the rotation transducer 
(potentiometer) and the voltage at the potentiometer 
rotor is fed into the strip chart recorder. When the 
control switch is in the "cal" position a dummy 
potentiometer (5 turn, 10,000 ohm, high resolution) is 
Switched into the circuit in place of the rotation trans- 
ducer. This facilitates the setting of calibration 
reference steps on the strip chart printout. 

Torque Transducer Signal Conditioner - The strain gage 
bridge of the torque transducer is connected to a Vishay 
Instruments, Inc., Model BA-4 Strain Gage Amplifier/ 
Signal Conditioner. The output from the signal conditioner 
is fed to the strip chart recorder. 

Motor Power Supply - The vane motor voltage is obtained 
from a Perkin Electronics Model TV RO40-15 D.C. power 
supply which supplies a regulated voltage variable from 
0-50 V.D.C. 

Strip Chart Recorder - Continuous torque rotation-time 
recordings is accomplished by a Honeywell Model 906 C 
"Visicorder" recording oscillograph. This machine has 


Six Channels available, two of which are used. Chart 
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speeds of .4, 2.0, 10.0 and 40.0 inches per second are 
available and the chart speed is controlled such that 
it is an accurate time scale. The speed was checked 


using a time mark generator. 





APPENDIX IV. - EQUIPMENT CALIBRATION 
AND CHECK QUT 
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A series of calibration and repeatability tests were conducted 


prior to core sample testing. 

Torque Measurement. - Linearity and repeatability of the torque 
transducer system were ascertained using the torque calibration 
apparatus described in Appendix III. The torque system was found 
to provide sensitivity of 0.001 in-1b with a.0.2 in-1b full scale 
setting and has a range of 0 to 10 in-1b available through gain 
steps in the signal conditioner. No measurable nonlinearity was 
found. 

Rotation Measurement. - The rotation measuring potentiometer was 
calibrated by mounting it on a large compass rose, with a pointer 
knob. 5 volts were put across the outside terminals and the 
voltage at the rotor terminal was measured as the knob was turned 
through 10° increments over the full throw of the potentiometer. 
The calibration factor was found to be 0.015 volts/degree or 0.859 
volts/radian, and the potentiometer was within 0.5% linearity. 
When the calibration factor is divided by the 1:367 vane shaft to 
potentiometer shaft rotation ratio the resulting voltage per vane 


shaft rotation ratio is 3.15] volts/radian. This calibration 


factor was used to obrain unadjusted vane rotation angles from the 


Strip chart read out. These angles were then adjusted by the 
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shaft twist at measured torque which was computed to be .0086 
radiangyin-Ib. 

Vane Motor Checks. - Motor speed versus input voltage was checked 
under loaded and no-load conditions. The loaded conditions were 
obtained by hanging weights from the calibration pulleys as in 

the calibration procedure. The motor speed was found load depen- 
dent at torques above its rated 10 in-oz. It was thus decided 
that only approximate (normal) rotation rates could be presetwith 
the motor input voltage. 

Time Scale Calibration. - It has been mentioned that the chart feed 
rate was used to obtain the times for the torque-rotation-time data. 
The accuracy of these rates for the recorder used was verified by 
connecting an electronic time work generator to a recorder channel. 
Check Out Tests. - A series of tests were conducted on three care- 
fully constructed remoided samples of submarine clays. Eachsample 
was aS nearly homogeneous as possible. The samples were at 
different water contents and their stiffness encompassed the 
expected range of stiffnesses of the hose hole samples. Each 
Sample waS teSted at 3 rotation rates with at least one rotation 
rate repeated twice. Each test was at an undisturbed location. 
These tests confirmed the capability of the apparatus to produce 
the required torque rotation-time-data over the entire rage of 
expected sediment shear strengths and at the wide range of speeds 


desired. In addition, the repeated runs (i.e. same sample-same 








speed) yielded the same torque versus rotation curves and thus 


further confirmed repeatability of the apparatus. 
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APPENDIX V. - X-RAYS OF SAMPLES TESTED 
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FIG. 42. - X-RAY PHOTOGRAPH, SAMPLE B-1, 19 FT 





FIG, 43, - X-RAY PHOTOGRAPH, SAMPLE B-1, 46 FT 
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FIG. 44. = X-RAY PHOTOGRAPH, SAMPLE B-1, 61 FT 
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FIG. 45. - X-RAY PHOTOGRAPH, SAMPLE B-1, 85.5 FT 
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FIG. 46, - X-RAY PHOTOGRAPH, SAMPLE B-1, 100 FT 





FIG. 47. - X-RAY PHOTOGRAPH, SAMPLE B-2, 22.5 FT 
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FIG. 48. - X-RAY PHOTOGRAPH, SAMPLE B-2, 40.5 FT 





FIG. 49. = X=RAY PHOTOGRAPH, SAMPLE B-2, 72.5 FT 
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FIG. 50. - X-RAY PHOTOGRAPH, SAMPLE B-2, 104.5 FT 





FIG. 51. - X-RAY PHOTOGRAPH, SAMPLE B-3, 55 FT 
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APPENDIX VI. - TORQUE VERSUS ROTATION 
AND LOG-LOG PLOTS OF G VERSUS TIME 
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